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Abstract 

In this paper we study a systematic and natural construction of canonical coor- 
dinates for the reduced space of a cotangent bundle with a free Lie group action. 
The canonical coordinates enable us to compute Poincare-Birkhoff normal forms of 
relative equilibria using standard algorithms. The case of simple mechanical systems 
with symmetries is studied in detail. As examples we compute Poincare-Birkhoff nor- 
mal forms for a Lagrangian equilateral triangle configuration of a three-body system 
with a Morse-type potential and the stretched-out configuration of a double spherical 
pendulum. 



AMS classification numbers: 70F07, 70G65, 53C80 



1 Introduction 

The theory of the reduction of Hamiltonian systems with symmetry is well developed 
[HE]- Although it is a classical subject and goes back to the pioneers of mechanics, a 
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modern theory was only established in the 1970's. The main idea can be summarized 
as follows. 

Let P be a symplectic manifold with a symmetry group G, and let J : P — > q* be 
an equivariant momentum mapping with respect to the coadjoint action of G on g*, 
where q* is the dual space of the Lie algebra g of G. Then, under some regularity 
conditions, the reduced space given by the quotient space P^ := J _1 (/i)/G A1 , where 
G^ is the isotropy group of £ q*, is a symplectic manifold. A G- invariant Hamilto- 
nian function on P can be reduced to a function on P^ which generates the reduced 
dynamics. 

The algebraic definition as a quotient makes it often difficult to explicitly con- 
struct the reduced space P^ and develop a good intuition for it. For example, P^ is 
not necessarily a linear space even if P is linear. This is the case, for example, in the 
ra-body problem: Although the translation reduced space is Euclidean, the reduced 
space of rotations is in general not linear [3] . But as the reduced space is a symplectic 
manifold it follows from the Darboux theorem that one can locally construct canon- 
ical coordinates so that the reduced space locally becomes a linear symplectic space. 
Such canonical coordinates are very useful. For example, they form the starting point 
of standard algorithms for the computation of the Poincare-Birkhoff normal form at 
an equilibrium point of a Hamiltonian system [H El El [8] . A Poincare-Birkhoff 
normal form is a main tool for the construction of center manifolds and the study 
bifurcations. As an example of the former application we mention the construction 
of the phase space structures which govern reaction dynamics induced by saddle type 
equilibrium points [91 [TO] . 

The main objective of this paper is the systematic construction of canonical coor- 
dinates for the reduced space P^ in the case where P is a cotangent bundle and the 
action of G on P is free, and illustrate how these coordinates can be used to compute 
Poincare-Birkhoff normal forms at the relative equilibria, i.e. the equilibria of the 
reduced system. This has numerous applications. To give one example we mention 
the construction of the phase space structures which govern the reactions in rotating 
molecules where the reaction dynamics is induced by saddle type relative equilibria 

p. 

We note that the computation of canonical coordinates for a reduced space of a 
symplectic manifold [12l[l3] in general or for specific cases such as, e.g., a cotangent 
bundle [144115] and more concretely for n-body systems [31 [TO] have at least implicitly 
been studied in the literature before. However, for obtaining the nonlinear terms of 
a Poincare-Birkhoff normal form these works have to be put into context, and a 
systematic study is missing. Also the work on the computations of Poincare-Birkhoff 
normal forms of symmetry reduced Hamiltonians is mainly restricted to Abelian Lie 
group actions. In this paper we present a systematic approach which covers both the 
Abelian and the non-Abelian case. 

In the following we give a brief review of existing literature related to this paper. 
In order to obtain canonical coordinates on the reduced space of a cotangent bundle 
with a Lie group action we follow the method given in [14] and [15] which take 
a Lagrangian respectively Poisson reduction point of view. A detailed survey on 
cotangent bundle reduction and its history can be found in [T7]. For the special 
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case of the three-body reduction, our main references are [18] and [3] to which we 
will come back in Sec. 01 As for the Poincare-Birkhoff normal form, one can find a 
detailed introduction in [19]. But for completeness, we give a brief explanation of 
the algorithm in Appendix [Cj One of the first applications of the normal form theory 
to reduced spaces of symplectic spaces with a continuous symmetry can be found in 
[2"U] where the symmetry group is the circle group. An application to the restricted 
three-body problem can be found in [21] and [22J, for instance. In [23] normal form 
computations are done at Lagrange points by using a splitting method. In another 
recent work [23] one can find a detailed study of normal form for planetary systems. 
Finally, a normal form at a relative equilibrium of a general dynamical system is 
given in [25]. A recent review of normal form theory in dynamical systems can be 
found in |S]. 

This paper is organized as follows. We start with a general review of the action 
of Lie groups on tangent and cotangent bundles in Sec. [2J This mainly serves to 
introduce some basic material and settle the notation. Sec. [3] comprises the main 
result of this paper which is a systematic construction of canonical coordinates for 
the reduced space of a cotangent bundle with a free action of a symmetry group. 
This includes the derivation of the reduced Hamiltonian in canonical coordinates, a 
detailed discussion of the case of simple mechanical systems (Sec. I3.3[) and of special 
cases like Abelian symmetry groups and systems with vanishing angular momenta 
in Sec. 13.41 and the Poincare-Birkhoff normal form of relative equilibria in Sec. 13.61 
Section [3] contains our first example which consists of the three-body reduction. We 
review in this section how to derive a reduced Hamiltonian in canonical coordinates 
in a way which does not depend on the choice of a body-fixed reference frame, i.e. in 
the language of Littlejohn and Reinsch [3] in a gauge independent way. In Sec. 14.31 
we consider a Lagrangian equilateral triangle relative equilibrium, we compute a 
Poincare-Birkhoff normal form at such configurations. The reconstruction of the full 
dynamics in the three-body case is addressed in Sec. 14.41 In Sec. [5] we study our 
second example which is the double spherical pendulum. After obtaining canonical 
coordinates for the reduced system, a normal form computation is done at the relative 
equilibrium given by the so called stretched out solution. Conclusions are given in 



2 Lie group actions on tangent and cotangent 
bundles 



In this section we recall the symplectic actions of Lie groups on tangent and cotangent 
bundles over a configuration space (mainly to introduce some notation). For the 
details, we refer to [261 E El [27] . 

Let G be a Lie group and let M be a manifold which is called the configuration 
space. Let the map 
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GxM->M 
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be a free action of G on M. We denote the left-translation which for a fixed g G G, 
maps s £ M to by L 9 . The derived maps of L 9 are denoted as follows. For s G M, 
(L g )* : T S M -> T gs M stands for the derivative map of L g , and (L s )* : T*M -> 
T*_ lg M stands for the pull-back map of L fl . Let g denote the Lie algebra of G. Then 
for £ G g, the corresponding infinitesimal generator or fundamental vector field 
at s £ M is defined by 



(^exp(tC) S ) • ( 3 ) 

i=0 



The G orbit through s G M is given by Gs = {55! j £ G} C M. The fundamental 
vector fields C,m are tangent to the orbits Gs for all s G M. Moreover, the tangent 
space T s (Gs) is spanned by the fundamental vector fields at s. 

If M = G, i.e. the action is the group operation of G, then the fundamental 
vector fields at g 6 G are given by 

C G (<7) = (^U, (4) 

where £ S g and i? g is the right-translation by g [I]. 

The coadjoint action of G on the dual space g* of its Lie algebra is defined as 

for g G G, n £ q* and ( G g. Here ( , ) stands for the pairing between a co-vector 
and vector and 

d 



AdgC 



dt 



( 5 (ex P (tC))5 _1 ) (6) 
t=o 



is the adjoint action. The action on M can be lifted to TM and T*M by the derived 
maps and both of the lifted actions are free when the action on M is free. The lifted 
action on the cotangent bundle is symplectic with respect to the natural symplectic 
structure on T*M and has a momentum mapping which is defined as follows: 
the momentum mapping J : T*M — > g* is given by 

(J(s, P ),Q = (p,(m(s)), (7) 

for all (s,p) G T*M. It is well-known that J is equivariant with respect to the action 
on T*M and the coadjoint action on g* [26] H]. 

If M is a Riemannian manifold with a Riemannian metric k which is invariant 
under the action of G, then the lifted action on TM is also symplectic with re- 
spect to the symplectic structure induced by the one on T*M and the corresponding 
momentum mapping L is defined by 

(L(s,v)X) =v T k( M (s), (8) 



for all {s,v) £TM|. 



3 Canonical coordinates in cotangent bundle 
reduction 



In this section we review the cotangent bundle reduction in the orbit reduction scheme 
of Marie [2S\ 127] with a coordinate-based approach. We then obtain canonical coor- 
dinates for the reduced space. A more detailed explanation of the notions used in 
the following subsection can be found for instance in [14|, |15j. 

3.1 Reduction of the equations of motion 

The shape space or internal space Q is defined as the quotient M/G. As we assume 
that the action of G on M is free it follows from standard theorems that M/G has a 
manifold structure and M — > Q is a fibre bundle [1] . Using the fibre bundle structure 
one can locally obtain a coordinate system on M by choosing a local coordinate 
system on Q x G. Let us assume that a point in Q x G has coordinates (q,g). 
The coordinates q are called shape coordinates or internal coordinates. Then by the 
decomposition [TJ] 

T(Q xG)^TQxGxg (9) 

a point in T*M has coordinates (q,q,g,g). Now consider the body angular velocity 
defined by 

e = (V0*& ( 10 ) 

where (L s -i)* denotes the differential of the left translation L g -i : G — > G, /i h> 
g h, at the unit element of G. The commonly used notion of body angular velocity 
comes from the fact that in the example where G is the Lie group 50(3) (see Sec. HJ) 
£ is indeed the angular velocity in a body fixed frame. Equation (fTQ|) is called the 
reconstruction equation as it can be used to find the full dynamics corresponding to 
the reduced one. We will comment on this in more detail in Sec. 14.41 for the case of 
three-body systems. 

One can see that the body angular velocity is invariant under the group action: 
for h £ G, define the curve m(t) = hg (t), then 

(L m -i)* rh = (I^-ih-i))* ((!*)* g) = £. (11) 

As £ is invariant under the group action, the coordinates (q, q, £) give a coordinate 
system on (TM) /G . 

If L : TM — > R is a regular Lagrangian function which is invariant under the 
action of G, then the function I : TM/G — > R given by 

i(q,4,0 ■= L (q,<i,g,g) (12) 

is well-defined. This is done by passing to the coordinates (<?, <?,£), and as the La- 
grangian L is invariant, it is possible to put L in the form of the function I in which 
the G coordinates disappear. Using I we can define momenta conjugate to q and £ 



6 



and 

•>=!• < i4 » 

respectively. Here r\ is called the body angular momentum, and by the chain rule, 

V = (L g )*p g , (15) 

where p g = dL/dg is the conjugate momentum of g G G Like body angular ve- 
locity the notion body angular momentum again comes from the context of reduction 
of rotational symmetries. 

Let H : T*M — > M be the Hamiltonian obtained from the Legendre transformation 
of the Lagrangian L |26J, i.e., in coordinates 

H(q, P q ,9,Pg) = qp q + gp g - L(q, q, g, g) . (16) 
As H is G invariant it induces a function h on T*M/G given by 

h(q,p q ,v) ■= H(q,p q ,g,p g ). (17) 
From the construction above one obtains for z = (q,p q , g,p g ) and £ G g, 

( J ( z ),() = ((.Pq,Pg),(M) 

= {(p q , P g),(0, (RgU)} 

= {Pg,(RgU) (18) 

= ({L g -l)*TI, (RgU) 

= <(MrO*??,C}. 

Here the first equality follows from the definition of the momentum map ([7]), the 
second equality makes use of the decomposition @, the fact that Qm is tangent to 
the group orbit which we identify with G and the third equality is clear, the 
fourth equality uses (|15|) and the final equality follows from the definition of the 
coadjoint action in ([5]) and ([6]). We thus obtain 

J(p) = (Ad g -i)* V (19) 

or equivalently 

r, = (Ad g yj(p). (20) 

Let stand for the coadjoint orbit through J(p) = fi 6 g* for some fixed \x G g*, 
i.e. 

O M = {(Ad p -0V| 5GG}cg*. (21) 

Then Eq. ([20]) gives that 

G O^. (22) 

Now consider the reduced space 

P^-J^iO^/G. (23) 
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By ([22]) . we conclude that if J(p) = /i for p = (q,p q , g,p g ) and some fixed [i £ q*, 
then for r\ = (Ad g )*J(p), we have (q,p q ,rj) G P^, i.e. (q,p q ,r]) are coordinates on 
the reduced space. The reduction to the space P^ we described is a coordinate- 
based form of the orbit reduction of Marie [27]. For n-body systems, the reduction 
procedure can be interpreted as passing to a body-fixed frame (cf. Seed]). In fact, the 
space P^ is symplectomorphic to the Marsden-W einstein reduced space J~ 1 (fi)/G^, 
where is the isotropy group of fi G q*. For us, the reduced space P^ and the 
coordinates (q,p q , rj) form the basis for defining canonical coordinates on the reduced 
phase space. 

3.2 Canonical coordinates 

The coordinates (q,p q , g,p g ) defined above are clearly canonical, whereas the coor- 
dinates {q,p q ,g,rj) as seen below are not canonical. By (fT3j) and (Q3|) we get the 
Poisson bracket equalities 

{qa,Va} = {Pg a ,r]a} = (24) 

on T*M for all a and a. Now recall the identification £ = (L g -i)* g, and let (ei, e{) 
be a basis of g, where I is the dimension of G. Choosing (q a ,qa, ga,t,a) in place of 
(qa,qa,ga,ga) as coordinates on TM requires one to obtain the dynamics in terms 
of the anholonomic frame (e±, ...,e;). We refer to [3] for a concise derivation of the 
Lagrangian and Hamiltonian in this context. We briefly give a derivation of the 
Poisson brackets in anholonomic frames in Appenix|XJ Using (|127p we get that 

{Va,Vb} = -labVc, (25) 

where 7^ fe are the structure constants given by [e a ,eb] = Jab e c- This is in fact the 
same as the (— ) Lie- Poisson bracket [2] on G (see also Appendix [B]) . As the Poisson 
structure on the coadjoint orbit is the reduced one from the Lie-Poisson structure 
on G, the discussion above suggests that if a canonical coordinate system (u, v) is 
chosen on such that 

rj = T](u,v), (26) 

then the coordinate system (q,p q ,u, v) becomes a canonical coordinate system on P^. 

For a more detailed discussion of canonical coordinates on coadjoint orbits, we 
refer to [30] among others. In our example of the three body-problem in Sec. H]the 
coadjoint orbits are the body-angular momentum spheres. In this case the canonical 
coordinates may be chosen as Deprit coordinates |31] on the body-angular momentum 
sphere (see ([85]) below). 

As for the dynamics, if one passes to the coordinates given in (|26|) . then the 
reduced Hamiltonian := h\p^ may be written in the form 

h^ = h^(q,p q ,u,v). (27) 

Then the equations of motion for the reduced system have the familiar form 

i = {zhf,} (28) 



for z = (q,p q ,u,v). 
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3.3 Simple mechanical systems 

The construction of canonical coordinates on the reduced phase space and the reduc- 
tion of the equations of motion can be described more explicitly in the case of simple 
mechanical systems. For a simple mechanical system, the Lagrangian is of the form 

L(s,s) = \s T ks-V{q), (29) 

where s is a coordinate system on M, A; is a Riemannian metric on M which is 
invariant under the action of G on M, and V is a potential function which is also 
invariant and hence, if (q, g) is a coordinate system on the local trivialization 

M = Q x G (30) 

then V depends only on the shape coordinates q. Let a : Q — > M be a local section 
of the fibre bundle, i.e. a is a right inverse of the projection n : M — )■ Q. If a point 
q G Q is given, then a point s G ^~ 1 (q) is of the form ga(q) with some g G G. Set 
r = a(q), then 

dr 

r = -q (31) 

by the chain rule. 

Let the body velocities be defined as 

v = (L g -i)J. (32) 

By using the Leibniz rule pQ and the definition of the fundamental vector field one 
obtains 

v = t[M(r) + f, (33) 

where £m is the fundamental vector field corresponding to £ = (L g -i)* g G q. The 
kinetic energy thus becomes 

K = -v T kv = -^(r) T 'k£ M (r)+ 6*(r) T kr + ^r T kr. (34) 
For £, rj G 0, set 

£ T I?7 = £M(r) T MM(r). (35) 

Then I is a left-invariant inner product on g [|. For G = 50(3) (respectively 
g = so(3)), I is the moment of inertia tensor (see Sec. H]). In order to decouple 
the kinetic energy in group and shape terms, the so called mechanical connection is 
introduced [2]. The mechanical connection A : TM — > g is defined by 

A(s,s) = I _1 L(s,s), (36) 

where L is the tangent momentum map given in ([8]), and I -1 : g* — > g is the linear 
map associated with the inner product I. At any point s G M, the tangent space to 
T S M may be decomposed into 

T S M = V S + H s , (37) 
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where V s is the tangent space to the orbit Gs, and H s is the space which is orthogonal 
to V s with respect to the metric k. A tangent vector w £ T S M may be written in 
this decomposition as 

w = ver s w + hor s w. (38) 

It turns out that [2] 

ver s w = [A(s,w)] M (s) (39) 

and 

L(s,hov s w) = 0. (40) 
With respect to this decomposition r may be written in the form 

r = ver r + (r — ver r). (41) 

If a new metric d is introduced by 

q T dq := (r — ver r) T k(f — ver r), (42) 

which is called the horizontal metric, and if we define the map Aq : TQ — > q by 

9 := ^4 r, (43) 

then after rearranging terms the kinetic energy assumes the form 

K=^(Z + A Q q) T I(Z + A Q q) + ^q T dq. (44) 

This compact form of the kinetic energy reflects the decomposition of it into vertical 
and horizontal energies. Finally the Lagrangian in coordinates (q, q, £) is given by 

Kq, Q, = \ (e + A Q q fl(Z + A Q q) + ±q T dq- V(q). (45) 
If the body angular momentum is defined by 

J = (Vi)»L, (46) 

then it is seen that 

J = I(Z + A Q q). (47) 
On the other hand, the conjugate momenta of q and £ are obtained as 

dl 

P q = ^=dq + A T Q 3, (48) 

and 

dl 

r ] = — = l^ + A Q q\ (49) 
respectively. Note here that n = J. Finally the Hamiltonian can be written as 

h(q,p q , J) = i J T r 1 J + (p, - 4; J) T d" 1 (p, - A T Q J) + y(g), (50) 

where d~ l denotes the metric on T*M corresponding to d. 
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3.4 Special cases 

(1) G is Abelian. Consider the case where the Lie group is Abelian, e.g. a torus 
group. The reduction strongly simplifies in this case. Since the coadjoint action is 
trivial the well-known identification [27] 

J- 1 (O fl )/G = T*Q (51) 

is obtained. This shows that one can take the coordinates (q a ,Pg) as the canonical 
coordinates on the reduced space, which is symplectomorphic to T*Q. An example of 
this situation is given by the translational motions of an n-body system [3] . We will 
explicitly illustrate the Abelian case for the example of a double spherical pendulum 
in Sec. [5j A detailed analysis of the Abelian case in the Lagrangian setting can be 
found in [T4"] . 

(2) Vanishing angular momentum. With the notation above, if ij = 0, then 
the coadjoint orbit is trivial as in the first special case. Then the reduced space 
is symplectomorphic to T*Q. A well studied example of a system with vanishing 
angular momentum is the so called falling cat problem |32| . 

(3) Generalized rigid bodies. Suppose that M = G. Then Q is just a point 

and 

J-\O ll )/G = O lt . (52) 

This occurs, e.g., for a rigid body, where the configuration space M is the rotation 
group SO (3) and the reduced space is the body-angular momentum sphere. 



3.5 Relative Equilibria 

A point on T*M is called a relative equilibrium point if its projection into the reduced 
space is a critical point of the reduced Hamiltonian. So we are interested in the 
equilibria of the function hn ■ Pn — > R for some fixed (i £ g*. We will give some 
criteria for relative equilibria for simple mechanical systems. 
From Eq. (|50p the equations of motion are obtained to be 

J = ~~ a d*dh/dJ J 

dh ,-i / at tn 
q= dp~ q = ( ^" A « J) ' (53) 

— = r^-A Q d-\p q -A T Q J) (54) 

V^=\j T I- 1 J + V{q) (55) 

is the effective potential. We note that for Abelian actions, the effective potential 
agrees with the so called amended potential [2j. The equations above can be deduced 



where 



and 
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from Appendix 1X1 and Appendix [Bj (See also [3l [33] and [M].) Then the conditions 
for having a relative equilibrium are 





= Aq J 


adj-i j J 


=0, 


oq 


=0. 



(56) 



In the examples of the three-body problem and the double spherical pendulum in 
Sections 0] and El the conditions will be given in a more explicit form. 

3.6 Poincare-Birkhoff normal form around a relative equi- 
librium point 

The Poincare-Birkhoff normal form is a main tool in dynamical systems theory. It 
allows one (if certain conditions are satisfied) to study the dynamics of a nonlinear 
system in the neighborhood of an equilibrium point by approximating it by a 'sim- 
pler' system. This has many applications, e.g., in the study of bifurcations and the 
computations of center manifolds [H El El E] ■ The simpler system is constructed order 
by order of the Taylor expansion of the original system at the equilibrium point by a 
suitable choice of coordinates at each order. For Hamiltonian systems, the coordinate 
transformation are sought to be symplectic. As the dynamics (i.e. the vector field) 
is generated by a Hamilton function the simplification can be described completely 
in terms of a simplification of the Hamilton function. There are well established al- 
gorithms which can be implemented on a computer and which allow one to compute 
normal forms to any desired order (see Appendix [Cj . As the starting point for these 
algorithms is a Hamiltonian system with canonical coordinates on the linear sym- 
plectic space R^ x R^ where / denotes the number of degrees of freedom it is crucial 
for the application of these algorithms to relative equilibria of symmetry reduced 
Hamiltonian systems to explicitly construct canonical coordinates on the reduced 
space as described in the subsections above. 

We will in this paper restrict ourselves to Poincare-Birkhoff normal forms at 
equilibrium points where the eigenvalues associated with the linearized Hamiltonian 
vector field ZD 2 H are purely imaginary. Here Z denotes the standard symplectic 
matrix and D 2 H is the Hessian of the Hamiltonian H. We will denote the eigenvalues 
by ±iwfc, k = 1, . . . , / . Assuming that the eigenvalues are independent over the field 
of rational numbers (i.e. in the absence of resonances), the Poincare-Birkhoff normal 
form yields a symplectic transformation to new (normal form) coordinates such that 
the transformed Hamiltonian function truncated at order no of its Taylor expansion 
assumes the form 

/ 

H NF (h, ...,//) = y^u k I k + h.o.t. , (57) 
fc=i 

where k = 1, . . . , / , are constants of motions which (when expressed in terms of 
the normal form coordinates) have the form 

h=pt + ql k = l,...,f, (58) 




Figure 1: Definition of the Jacobi vectors s x and s 2 and the corresponding angle <fi. 



and -ffNF is a polynomial of order no/2 in 1^, k = 1, . . . ,/ and hence of order no 
in p and q (note that only even orders no of a normal form make sense). The 
algorithm to compute this transformation is sketched in Appendix O We will apply 
it to the examples of relative equilibria of a three-body system and a double spherical 
pendulum in Sec. 14.31 and Sec. 15.31 respectively. 



4 Three-body systems 

In this section we review the reduction of a three-body system for which we then 
write the reduced Hamiltonian in canonical coordinates following Sec. [3l We also 
comment on the reconstruction of the full (unreduced) dynamics. As an example we 
discuss a triatomic molecule with a Morse-type potential for which we compute the 
Poincare-Birkhoff normal form about an equilibrium point given by an equilateral 
triangle configuration. 

4.1 Reduced equations of motion 

Consider a system of three bodies with masses mi,m 2 ,m 3 and position vectors 
xi,xi,xi G M 3 , respectively, without external forces acting on the three bodies. 
The symmetry of overall translations can be reduced by introducing mass-weighted 
Jacobi vectors which are defined according to 

si = V/^(xi-x 3 ), 

— mixi + m 3 x 3 

S2 = V/M X 2 1 ), 

mi + m 3 

where 

mim 3 m 2 (mi+m 3 ) 

/*i = : . = : . (59) 

mi + m 3 mi + m2 4- m 3 

are reduced masses (see Figure [1]) . 

Excluding collinear (and hence also collisional) configurations we obtain the six- 
dimensional translation-reduced configuration space 

M = {s = (si,s 2 ) : Asi + us 2 ^ for all (A, n) £ M 2 \{0}} C I 3 xR 3 . (60) 



Proper rotations g G SO (3) act on M in the natural way 

s(si,s 2 ) = (gsi,gs 2 ). (61) 

On M this action is free and it thus follows that the shape space 

Q := M/SO(3) (62) 

has a manifold structure which turns out to be diffeomorphic to M.^_ = {(x,y,z) G 
M 3 : z > 0} [12] • The canonical projection M — > Q defines a principal bundle with 
structure group SO (3). This principal bundle is trivial [18] and has holonomy group 
SO (2) yielding a geometric reduction [35J. 

The Lie algebra g = so (3) of G = SO (3) can be identified with M 3 where the Lie 
algebra structure becomes the vector product 'x'. By using the bi-invariant inner 
product on g, or equivalently the dot product on M 3 , one can identify g* also with 
M 3 . With these identifications the fundamental vector field corresponding to ( G g 
at s = (si, S2) is 

C(s) = (Cx Sl ,Cxs 2 ). (63) 
The momentum mapping L : TM -4 g*, following (|8|), is given by 

L = si XS1+S2 x s 2 . (64) 

After choosing a body-fixed frame one can obtain the corresponding body-fixed 
Jacobi vectors by 

Si=gr i: £ = 1,2, (65) 

where g G SO(3) is the matrix relating the body- fixed frame and the space-fixed 
frame. As g depends on three coordinates, e.g. Euler angles, there are three shape 
space coordinates q a , a = 1,2,3, remaining to parametrize the two vectors ri and 

The kinetic energy is given by 

K = \jZm^ =\Y,si ■ (66) 

i=i i=l 

The corresponding metric k thus is Euclidean. Defining body velocities according to 

Vi=g T Si, £ = 1,2, (67) 

(cf. (|32p ). and using the shape coordinates and their time derivatives one can rewrite 
the body velocities as 
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The body angular velocity £ (see (JTOj) ) is the vector in R 3 corresponding to S G Q 
given by 

= = 5 9, 

which is the reconstruction equation (|10p . Then one has 

3 



Q = l 



<9r 



"9a 



(69) 



which corresponds to the general expression (|33p . Since the moment of inertia tensor 
I is given by 

lu = ri x (u x n) + r 2 x (u x r 2 ), (70) 

Al 

for ueK 3 , the mechanical connection Aq 



A 2 
A 3 



is obtained to be 



(n x — ) + 1 r 2 x — . 

dq a dq a 



Then the kinetic energy becomes 



3 3 



where 



a=l 



a,B=l 



dr{ T dri 
^ dq a dqp ' 



Using that the horizontal metric is 

da/3 = h a p — I Ap 
(see (|42p ) allows one to write the kinetic energy in the compact form 

1 3 1 

a,P=l 

Following (|46|) the body angular momentum is given by 

J = g L = ri x vi + r 2 x v 2 . 
Then by Eq. (j69|) and one has 



Q = l 



The conjugate momenta are given by 

dK 



9e 



I-(£ + 5^A a g Q ) = J, 



(71) 



(72) 



(73) 



(74) 



(75) 



(76) 



(77) 



(78) 



a=l 
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and 



Pa = = ^d a pqp + 3 T A a . (79) 
° Qa 0=1 

Thus the Hamiltonian takes the form 
1 1 3 

h=-J T r l J + - ^ d af3 ( Pa -J T A a )(pp-J T A ) + V, (80) 

<x,P=l 

where V = V(q±, q2, (73) is the potential. 

Let us now make give explicit expression by introducing coordinates. As the shape 
coordinates (<Zi, <72j 93) we choose Jacobi coordinates (ri,r2,(f>) which are defined as 
® 



V r i • r i> 



T2 



V r 2 • *2, <fi = cos 1 (r 1 ■ r 2 /{r 1 r 2 )), < cf> < vr, (81) 



(see Fig. [[]). Choosing then the axes Xb,Vb, %b of a body-fixed frame according to the 
so called xxy-gaug^i] shown in Fig. [21 one obtains for the moment of inertia tensor, 
metric and mechanical connection [3] 



r\ sin 2 (j) 
-r\ sin 4> cos < 




— r\ sin (j> cos < 
000 

r{ + COS (f 








(82) 



[d 



1 
1 





and 



(0,0,0), 



(0,0 



' r\ + r\ 



(83) 



(84) 



respectively. 

Putting the results above together the Hamiltonian in terms of Jacobi coordinates 
becomes 



h{ri,r 2 ,<f>,Pi,P2,P3, J) 



1 r r\ + 7*2 COS 2 (ft j2 _|_ 2 COS (j 

2 r 2 r| sin 2 (j) 1 r 2 sin < 

„2 



Ji J2 H — 2^2 "h 



r 2 + r| 3 



+ Vl+Pl+ 1 2 2 2 CP3 



r 2 + r\ 



J3) 2 } + V( ril r 2 , 



where J = (Ji, J2, J3) . Here ||J|| is conserved so the coadjoint orbit is the body 
angular momentum sphere 5 2 (|| J||). One choice of canonical coordinates on 5 2 (||J||) 
are the so called Deprit coordinates which are defined as [31] 



(Ji, J2, J3) = (v, V V 2 — v 2 sinu, 



\/ r 2 — v 2 cos u) 



(85) 



1 We note that the choice of a body-fixed frame corresponds to the choice of the local section a of the 
fibre bundle M — > Q in Sec. 13.31 The gauge theoretical interpretation of this choice is studied in great 
detail in [3J. 
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Figure 2: Definition of a body-fixed frame according to the xxy-g&uge. 

are chosen, where r = ||J|| (see Fig. [3]), then the reduced Hamiltonian h\gar r \ becomes 

1 rf + r% cos 2 (/> 2 2 cos /-= „- . 

h r {n,r2, <p,Pi,P2,P3, u,v) = -{ — . 2 u +i— rWr "» smu 

2 rf r| sin <p rf sm 

1/2 2\ • 2 , 1 ( 2 2\ 2 ,2,2 

+ — (r — u J sin u+~2 2V ~ v ) cos n +Pi+y , 2 

r l r l + r 2 



2 2 

+ ri ^r 2 (p 3 - / 2 W r 2 -v 2 cosu) 2 } + V. 

r l r 2 r l + r 2 



4.2 Lagrangian equilateral triangle configurations 

We now consider a three-body system with a Morse-type potential given by 

3 

V= exp(-2(r ii -d ))-2exp(-(r ii -do)), (86) 

1<«<J<3 

where r,j is the distance between the ith and the jth particle. The parameter do 
determines the side length of the equilateral triangle at which the potential has a 
minimum. In Jacobi coordinates one has 



ria = (87) 



^23 



ri2 
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Figure 3: Coordinate lines on the angular momentum sphere of (u, v) defined according to 

am. 



A Lagrangian equilateral triangle relative equilibrium is a planar motion where the 
shape is a constant equilateral triangle. The angular momentum is orthogonal to the 
plane of the motion. The mass-weighted Jacobi vectors are then of the form 




(90) 



where the parameter b is determined by the magnitude of the angular momentum r or 
conversely, choosing a value for b determines r. The corresponding Jacobi coordinates 
(rf,r|,(/> e ) are easily computed using (|8ip. 

Now we find the values of the other coordinates and the parameter r at the 
equilibria specified by b. Following ([80]) a relative equilibrium satisfies [33] 



J x (IT 1 - J) = 0, (91) 
Pa = J-A Q , (92) 

±. ( \^ 3 + V) = 0. (93) 

By Eq. (|9ip J is an eigenvector of I -1 at a relative equilibrium point. For a Lagrangian 
equilateral triangle relative equilibrium we know that in the xxy-gauge J is pointing 
in the z-direction of the body frame. Hence J = (0,0, r). From (|85p we find the 
corresponding canonical coordinates (u e ,v e ) = (0,0). Inserting J = (0,0, r) and 
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Figure 4: (a) The magnitude r of the angular momentum as a function of the size of the 
Lagrangian equilateral triangle parametrized by b (see ( 190]) ). (b) The energy as given by 
the effective potential V e g = \ Z> i + V (see (BID ) f° r the Lagrangian relative equilibria as 

z r 1 +r 2 

a function of r. 



using the block structure of the inertia tensor in ([82]) Eq. ([93]) reduces to 



9,1 r 2 



— (o^— -i+V) = 0. (94) 



We use this equation to find the magnitude of the angular momentum r for a La- 
grangian equilateral triangle (rf,rf,0 e ) specified by a given parameter b in ([90]) . 
Finally, inserting (|84|) in (|92p the conjugate momenta are obtained to be 



rf=Pi=0, ri= ^ ) if ( ^ r. (95) 



Figure Hti shows the magnitude of the angular momentum r as a function of 
the Lagrangian equilateral triangle specified by b. One sees that for a given value 
of r, there are two (or no) Lagrangian equilateral triangle of different size. The 
corresponding energies given by the effective potential V e g = h J, $ + V at these 

r l+ r 2 

equilibria are shown in Fig. Hb. For a given value of r, the smaller Lagrangian 
triangle has the smaller energy. 

As mentioned above at a relative equilibrium point the body angular momentum 
vector J is an eigenvector of I -1 . So, when looking for relative equilibria in gen- 
eral one would like to diagonalize I -1 which is possible when passing to a principal 
axes frame. We note that the corresponding shape coordinates are called Draght's 
coordinates [21 [T5]. As we were only interested in Lagrangian equilateral triangle 
configurations in this paper, the commonly used Jacobi coordinates were also useful 
in the study of these relative equilibria since I -1 is diagonal in the third compo- 
nent which corresponds to the direction of the body fixed angular momentum in the 
a^y-gauge. 
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4.3 Normal form around Lagrangian equilateral triangle 
relative equilibria 

We now apply the procedure explained in Appendix [C] to compute the Poincare- 
Birkhoff normal form around the Lagrangian equilibria. We choose unit masses, the 
parameter do in the Morse potential in (j86[) equal to 6 and the parameter b specifying 
the side length of the Lagrangian equilateral triangle in (|90p equal to 6.5. This gives 
the Jacobi coordinates 

e e , u 6 - 5 

Qi = n = v^io = ^= > 

Solving dMD for r we find r e = 19.8302179854. 

The momenta conjugate to the Jacobi coordinates are 

Pi = 0, 

Pl = 0, 



(97) 

-r = Q Q1fiin8QQ9.7 

{r{) 2 + (rf) 5 

The eigenvalues of the matrix $D 2 h r that gives the linearized vector field are 



^ = WT7^2 r = 9 - 9151089927 - 



±iwi = ±i 0.2362174000, 
±iw 2 = dbi 0.4693542718, 

(98) 

±iw 3 = ±i 1.1749259437, 
±iw 4 = dbi 1.1984363284. 

So we can immediately read off that the equilibrium is of elliptic linear stability. 
We note that a well established method for determining the stability of reduced 
systems is the reduced energy-momentum method which was introduced in [36]. For an 
application to the three-body problem, see also |37j . The reduced energy-momentum 
method does however not provide a means to compute higher order normal forms as 
we will do now following AppdenixlO 

Since we are only interested in demonstrating the basic principle of a normal form 
computation we will restrict ourselves to the normal of order 4. We start from the 
fourth order Taylor expansion of the Hamiltonian h r at the relative equilibrium. It 
has 212 nonvanishing terms and we refrain from writing them down. The symplectic 
matrix M which yields the linear symplectic transformation (|142|) after which the 
quadratic part of the Hamiltonian assumes the form 

4 
k=l 
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(see (|144p ) can be denned as 

M = [ciRe vi, C2Re V2, C3Re V3, C4Re V4, cxlm vi, C2lm V2, C3lm V3, c^lm V4] (100) 



where the column vectors are the real and imaginary parts of eigenvectors v& of 
$D 2 h r for the eigenvalues iw^ with coefficients 



Ck 



VRe v k ■ Jim v k 



1,...,4. 



(101) 



We find M 



fo 



V 





5952500442 
0.5952500442 





3.2144619462 



0.2245619939 












0.5952500442 
-0.5952500442 





0.6462635772 





-0.6459181965 
-0.6459181965 

















0.2590186725 



0.1406084178 
-0.1406084178 










4.4531132913 










-0.2590186725 


0.6993747200 
-0.6993747200 




(102) 

Following the next step in the procedure described in Appendix O we find that the 
normal form of order 4 is given by 



\ 









-0.7740918318 
-0.7740918318 
/ 



J*) 



NF =2.1181531267 + 0.2362174000 ii + 0.4693542718 I 2 + 1.1749259437 I 3 + 1.1984363284 / 4 

- 1.4978871558 if - 7.7221894156 Iil 4 - 0.9580186364 l\ + 6.4183166825 Ii_I 3 - 8.1361397396 I3I4 (103) 

- 0.8641444715 l| - 0.2175152611 - 0.2069751620 I2I4. - 0.1815241432 I 2 I 3 + 0.0089977794 l\ . 



The Hamiltonian is obtained from the general approach described in this paper. 
It yields an integrable nonlinear approximation of the 3-body problem reduced by the 
non-Abelian symmetry group SO (3) which can be used to study the motion in the 
neighborhood of the Lagrangian equilateral relative equilibria. Higher order terms 
can be obtained following the procedure in Appendix ICl 



4.4 Reconstruction of dynamics 

Generally speaking, when a curve in the reduced space is given finding the actual 
curves in the full space is the problem of reconstruction. This topic is well developed 
and we refer to |38j for the details. We briefly sketch how the full dynamics can be 
computed in the three-body reduction case. 

Consider a curve c M (t) = (ri(t), r 2 (t), <t>(t),Pl(t),P2(t),P3(t), u(t), v(t)) in the re- 
duced space. Then what is the corresponding curve in the full space of which the 
projection is c M ? Firstly, we can find J = (Ji,J2,Js) by (|85p. Then by Eq. (|77|) it 
is easy to obtain the angular velocity vector £ or the corresponding matrix S. After 
that one has to solve the reconstruction equation (|68p which is not a trivial task 
because of the non-Abelian structure of SO(3) (see the explanations given in |38|). 
A solution g gives the SO (3) coordinates, and finally by Eq. (fT5|) their conjugate 
momenta are obtained. For a more detailed discussion of the reconstruction in the 
case of the three-body problem, we refer to [39] , 



Figure 5: The double spherical pendulum. 



5 The double spherical pendulum 

In this section we study the reduction of the double spherical pendulum (for a more 
detailed survey, we refer to [10]). We again introduce canonical coordinates on the 
reduced space in the light of Sec. [3l We use these to compute the Poincare-Birkhoff 
normal form at the relative equilibrium given by a so called stretched-out solution. 

5.1 Reduced equations of motion 

Consider two coupled spherical pendula with masses mi and m.2 and position vectors 
si and S2 defined as in Fig. [5j moving without friction under the influence of a 
gravitational force — ak where a is a positive constant and k is the unit vector in the 
z-direction. 

If the lengths of si and S2 are l\ and 1%, respectively, then the configuration space 
is M = S 2 (h) x S 2 (l 2 ). The Lagrangian is 



The system is invariant under rotations around the z-axis. So the symmetry group 
is the Abelian group S 1 whose action on M is given by 



L(si,s 2 ,si,s 2 ) 




(104) 
(105) 




(106) 
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where gg is the rotation by the angle about the z-axis. We can identify the Lie 
algebra of S* 1 with span(k). An element of the Lie algebra is then an angular velocity 
vector of the form wk with w£R and the corresponding fundamental vector field is 
oj (k X Si, k X S2). For the angular momentum, we find according to ([8]) 

(L(si,s 2 ,si,s 2 ),o;k) = w(misf (k x si) + m 2 (s 1 + s 2 ) T (kxsi + kx s 2 )) (107) 

or 

L = k T (mi(si x si) + m 2 (si + s 2 ) x (si + s 2 )) k . (108) 

If the body frame is chosen such that the x-axis coincides with and we introduce 
polar coordinates (r, 0) in the xy-plane then we obtain for the body-fixed position 
vectors 



n = (nA-^lf-rf), (109) 

r 2 = (r 2 cos (p, r 2 sin 99, - yft% - r|) . (110) 

Note that through the choice of the sign of the square roots in the last components 
these equations are restricted to downward pointing configurations. As rx,r 2 ,<p are 
invariant under the group action we can take them as shape coordinates on the 
three-dimensional shape space Q = S 2 {1\) x S 2 (l 2 )/S 1 . 
The moment of inertia tensor 

I = mi ||rf- 1| 2 + m 2 ||(n + r 2 )^|| 2 , (111) 

where rj~ is the projection of ri onto the xy-plane, can be written in terms of the 
shape coordinates as 

I = (mi + m 2 )r\ + 2m 2 rir 2 cos (p + m 2 r 2 . (H2) 
Accordingly, we get for the mechanical connection 

m\m 2 r 2 sinip mim 2 rismip mim 2 r 2 (ri cos tp + r 2 ) 

A-rx = J 1 A-2 = J j Ay = 1 , (113) 

and the entries of the matrix d which gives the horizontal metric are 
Z 2 (mi+m 2 ) m\m\r\ sin 2 ip 



dt 



— rf) (mi + ra 2 )r\ + 2m 2 rir 2 cos ip + m 2 r\ ' 



if I 1 2m 2 m 2 sin ra 

112 = 7rm 2 cos ip + r x r 2 {—===—=== + 



2 \ \ - r\^Jl\ - r\ (mi + m 2 )rf + 2m 2 rir 2 cos 99 + m 2 r| 

, 1 / 1 2m?m 2 r 2 (ri cos 99 + r 2 ) \ . 

113 = ^m 2 r 2 -1 + 7 ; v 2 : ; 2 Sinyj, 

2 \ (mi + m 2 )r^ + 2m 2 rir 2 cos 99 + m 2 r2/ 



,o 9 9-2 

/ <o mf m 2 rf sm <z> 

^22 = m 2 1 " 



2Z| — 2r 2 (mi + m 2 )r 2 + 2m 2 rir 2 cos ip + m 2 r| 

m\m\r\r 2 {r\ cosy + r 2 sin 99) 

(mi + m 2 )r 2 + 2m 2 rir 2 cos 99 + m 2 r 2 ' 

1 9/ 2m 2 m 2 (r\ cos 0? + r| 

^33 = 77^2^2 11 



2 y (mi + m 2 )rf + 2m 2 r\r 2 cos y? + m 2 r^ 
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The left-action of the group on the tangent bundle is trivial, so J = L. The 
conjugate momenta of the shape coordinates are given by 

dL 3 

Pa = -K— = ^Jftj^/S + J A a . (114) 
° Qa 0=1 

Thus the Hamiltonian takes the form 

3 

T " 2 



1 1 3 

h = -r 1 J 2 + - d a P( Pa -JA a )(pp-JAp) + V, (115) 



a,B=l 



where V = —mi a ytf — r\ — mi a — r\ + \Jl\ — r%) is the potential. Observe 
here that as L is conserved, J is also conserved and can be viewed as a parameter. 
The reduced equations of motion are 

dh d 1 (116) 

Pa = ~ P = -^(-{r 1 J 2 + d°*(p a - J Aa) {pp - J Ap)} + V(q)) , 
uq a oq a z, 

where we denote by r the z-component of the conserved angular momentum J. 



5.2 Relative equilibria 

Because of the triviality of the coadjoint action the conditions to have a relative 
equilibrium (j56|) reduce to 

Pa = J A a , 

As shown in [40J there are two types of relative equilibria: the so called cowboy 
branch and the stretched-out solution that we will concentrate on in the following 
and which is shown in Fig. [6j For a stretched-out relative equilibrium, we have (p = 0. 
We find ri and ri from solving (|117p (using the computer algebra program Maple). 
The corresponding momenta are obtained from (jl 14)) . We will in the following choose 
all parameters to have unit values, i.e. m\ = l\ = mi = li = a = 1. The energy of 
the stretched-out relative equilibrium as a function of J for this choice of parameters 
is shown in Fig. [7k. 



5.3 Normal form around stretched-out relative equilib- 
ria 

The stretched-out relative equilibria are known to be stable [3D]. In agreement with 
this result we find that the eigenvalues of the matrix 2D 2 h r associated with the 
linearized vector field of the reduced system at the stretched-out relative equilibria 
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Figure 6: A stretched-out relative equilibrium solution of the double spherical pendulum 
in which the two masses are aligned with the point of suspension and move along circles. 




r 



Figure 7: (a) Energy of the stretched-out relative equilibrium as given by the effective or 
amended potential V e s = J 2 + V(q) (see (11171) ) as a function of the angular momentum 
r. (b) Frequencies of the stretched-out relative equilibrium as a function of r. 
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are purely imaginary. The frequencies u^, k = 1,2,3, are shown in Fig. [7)3 as a 
function of r. 

For the normal form computation, we consider the relative equilibrium point 
which has r = 1. For the position of this relative equilibrium, we find 

q\ = r \ = 0.4425598655 , 
q\ = r e 2 = 0.5656579210 
q % = m e = , 

3 I (H8) 
Pi = o, 

P2 = 0, 

p 3 = 0.4704091824 . 

The Taylor expansion of the reduced Hamiltonian to order 4 has 186 nonvanishing 
terms at this relative equilbrium. The eigenvalues associated with the linearized 
vector field are 

±iu>i = ±i 1.2572610531 

±iw2 = ±i 1.4864684140 (119) 
±\oj 3 = ±i 2.6603546311 

We define the symplectic matrix M which transforms the quadratic part of the Hamil- 
tonian to the form Ylk=i + Ik) analogously to (|100|) in Sec. 14.31 We find 



M 



0.3720476175 -0.3116300067 0.4202281883 

-0.5712604029 -0.3238407419 -0.5369033436 

2.2390869882 -0.2974563002 -2.2029538091 

-0.8426837965 -1.6576611267 1.8964495311 I 

0.3873035344 -1.4365053634 -1.4081719163 / 

-0.6854444155 0.0612664420 -0.2510237658 / 



(120) 

Following the procedure described in Appendix [C] we get for the 4th order normal 
form 

h^ F = - 2.2056999577 + 1.2572610531Ji + 1.4864684140/ 2 + 2.6603546311/ 3 

+ 0.0467015469i 1 2 - 5.8213832524/i/ 3 + 0.0875786340/f + 0.1772800788/i/ 2 
- 0.0932948515/ 2 h - 0.0419637147 if . 

(121) 

The Hamiltonian yields an integrable nonlinear approximation of the double 

spherical pendulum reduced by the Abelian symmetry group SO (2) in the neigh- 
borhood of the stretched-out relative equilibrium. Similarly to the 3-body case in 
Sec. l4.3l the normal form of the reduced system is obtained from the general approach 
in this paper which demonstrates the effectiveness and generality of the approach. 



6 Conclusions and Outlook 



In this paper we provided a general perspective on the construction of canonical co- 
ordinates for the reduced spaces of Hamiltonian systems given by cotangent bundles 
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with a free Lie group action. The general approach presented in this paper allows 
one to treat the reduction of Abelian and non-Abelian group actions on the same 
footing. The case of simple mechanical systems was studied in detail. The approach 
was illustrated for a 3-body problem and the double spherical pendulum which in- 
volve non-Abelian and Abelian symmetries, respectively. We used the canonical 
coordinates to compute the Poincare-Birkhoff normal forms at the relative equilibria 
given by the Lagrangian equilateral triangle configuration in the 3-body problem and 
the stretched-out solution of the double spherical pendulum. The Poincare-Birkhoff 
normal form gives a nonlinear approximation of the local dynamics of the reduced 
system in the neighborhood of the relative equilibria. This goes beyond the well 
established reduced energy-momentum method which only give the linear stability 
of the relative equilibria [36]) and enables one, e.g., to give nonlinear approximations 
of the center manifolds of relative equilibria. The use of a Poincare-Birkhoff normal 
form for the computation of the center manifolds of saddle type equilibria has in 
recent years been demonstrated in the study of reaction type dynamics [10]. The 
study of this paper allows one to carry over these results to the case of saddle type 
relative equilbria which induce reaction type dynamics in rotating molecules In 
this context also the reconstruction of the full dynamics from the reduced one which 
we illustrated for the example of the 3-body problem is of importance. 

In this paper we excluded Lie group actions with isotropy which have been studied, 
e.g., in [13] or [33J and [16J for the case of n-body systems. Our future studies concern 
how isotropy can be incorporated in the approach presented in this paper to do, e.g, 
a normal form analysis for 3-body systems with linear equilibrium configurations. 
Another related problem is the development of a normal form algorithm which is 
coordinate independent or in the jargon of [3] gauge independent. This problem is 
considered in [41] and its generalizations to cotangent bundles in general seems worth 
studying. 

Recently considerable progress was made in nonholonomic mechanics (see, e.g., 
[291 |4"2] and the references therein). It might be possible to use these techniques 
to develop a non-canonical Poincare-Birkhoff normal form related the results of this 
paper. 

Finally we mention that it would be interesting to transfer the results of this 
paper to quantum mechanical systems. Analogously to the Poincare-Birkhoff normal 
form of an equilibrium point there is a quantum normal form built on the symbol 
calculus of pseudo differential operators by which one can locally approximate a 
quantum Hamilton operator. In the case of elliptic equilibria this allows one to 
compute quantum energy spectra with high precision. For saddle type equilibria, the 
quantum normal form can be used to compute efficiently quantum reaction rates and 
the associated Gamov-Siegert resonances (see [43 \ |4"4"] and also the references therein 
for quantum normal forms in general). It would be interesting to transfer these 
results to relative equilibria of rotational symmetry reduced molecular systems. The 
dependence of the quantization of quantum reaction rates in the hydrogen exchange 
reaction as a function of the angular momentum has, e.g., been studied in [45] using 
ab initio quantum computations. It would be interesting to compare these results 
to a quantum normal computation. The geometric approach for quantum 3-body 
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problems presented in [18] could be very useful for this purpose. 

A Poisson brackets in an anholonomic frame 

We give the Poisson brackets in an anholonomic frame, the details can be found in 
[3]. Let M be a manifold and s\, s n be a local coordinate system on M. Consider 
a regular Lagrangian L : TM — > R which then is a function of the coordinates 
s\, s n , s\, s n , and a local frame X\, X n on M with 

n 



where di = d/ds{. From 



we obtain that 



S = ^2sidi = ^2 V i X i' ( 123 ) 
3=1 j=l 



dsj 
dvj 



a). (124) 



Let L(si,Vi) = L{si,Si). Then the conjugate momenta Pi = §jr; and = are 
related by 

7T, = a{pj . (125) 

If the Poisson bracket of two functions / and g on M in terms of the coordinates 
si, s n ,pi, ...,p n has the canonical form 

~{ d s i dPi dPi ® s i ' 

then using Equations (|124p and (|125[) the Poisson brackets of the coordinates 
si, ...,s n ,7Ti, ...,n n are given by 

n 

{Si,Sj} = 0, {Si,TTj} = aj, {xj,^} = -^CyTTfc, (127) 

k=l 

where [X^X,-] = J2k=i c ij^k- Then the Poisson bracket f|126j) becomes 

B Lie-Poisson structures 



We recall some basics of Lie Poisson structures on Lie groups. For the details we 
refer to [2]. Let G be a Lie group and g be its Lie algebra with the Lie bracket 
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[ j ] : § x g -> g. The dual space g* of g is a Poisson manifold with either of the 
following brackets 

</. *> ± 00 -*(/.. [5.5;] )■ d») 

Here g** is identified with g in the sense that 5f /5/j, G g is defined by (f , Sf/S/j) = 
D/(/i), where D is the derivative. Let B = e\, ...,e m be a basis of g and e\, —,e^ 
be the dual basis of B in g* with the corresponding coordinates £i,...,£ TO 011 8 an d 
/ii, ...,^ m on g*, respectively. Then the Lie Poisson bracket is given by 

{f,k} ± = ±f, a ^ c 6 / 6 / , (130) 

where [e ffl , e b ] = 7^ e c . 

Let H : q* — > R be a Hamiltoninan function. Then the general Lie Poisson 
equations determined by F = {F, ff} read 

A = ±ad «f/5/^> ( 131 ) 
where ad^ : g — > g is the linear map ( — > [£ , (], and ad| : g* — > g* is its dual. 

C The Poincare-Birkhoff normal form 

Let H2 denote the quadratic Hamiltonian which gives the linearized Hamiltonian 
vector field at the equilibrium. One says that a Hamiltonian H is in normal form if 
H Poisson commutes with its quadratic part, i.e. 

m m ._4- ( 9H 2 dH dHdH 2 

- X, i WkWk ~ Wk ^) " 0, (132) 

where / is the number of degrees of freedom. In general H is not in normal form. 
However, for any given order no of the Taylor expansion of H one can find a symplec- 
tic transformation to new phase space coordinates in terms of which the transformed 
H truncated at order no is in normal form. This symplectic transformation is con- 
structed from a sequence of the form 

(ft, Pi) = z = z(°) 1 — y z« 1 — y 1 — y 1 — y ... 1 — y z( n °) , (133) 
where z^ is obtained from z*-™ -1 ) by means of a symplectic transformation 

z (n-l) ^ z (n) = ^ z («-l) (1 34 ) 

generated by a polynomial W n (z) of order n, i.e. 

W n G W n := span {g^ 1 . . . qj'pf 1 . . .pj' : |a| + |/3| = n} . (135) 

Here \a\ = ^{ =1 «fc, \P\ = Ylk=i@k- More precisely, the <f>w n in (|134p denote the 
time-one maps of the flows generated by the Hamiltonian vector fields corresponding 
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to the polynomials W n (see [H] for the details). The maximum order no in f)133j) is the 
desired order of accuracy at which the expansion will be terminated and truncated. 

Expressing the Hamiltonian H in the coordinates z^ n \ n = 1, . . . , no, we get a 
sequence of Hamiltonians , 

H = H<® -> -> i^( 3 ) . . . -> fl< n °) , (136) 

where for n = 1, . . . , n , H^ n \z^) = fft"- 1 )^"- 1 )) = ^("^'(^zW), i.e. 

(») = H (n-1) Q ^-1 (137) 

To avoid a proliferation of notation we will in the following neglect the superscripts 
(n) for the phase space coordinates. 

In the first transformation in (|133|) we shift the equilibrium point zq to the origin, 
i.e. z i y 4>w 1 (z) := z — zq. This gives 

H«(z) = ff(°)(z + z ). (138) 

The next steps of the normal form procedure rely on the power series expansions of 
Hi"), 

oo 

H^(z)=E + Y,Hi n \z), (139) 

where the Hg are homogenous polynomials in W n : 

Lrfn)^ \ ai,...,Q!f,/3i,...,/3f Q ai A /?/ , 1 

//; ^ ai \...a f W...p f \ q i ■■■ q f P " - P f (M0) 

|q| + |/3|=s j j 

For n = 1, the coefficients in (|140p are given by the Taylor expansion of at the 
origin 



H {1) - TT H m (£ 



(141) 

z=0 

For n > 3, the coefficients in (|14(jp are obtained recursively. For n = 2, i.e. the 
second step in the sequence of transformations (|133p . the coefficients in ()140p are 
determined by a linear transformation of the phase space coordinates according to 

z h-> <pw 2 {z) ■= Mz. (142) 

Here, M is a symplectic 2/ x 2/ matrix which is chosen in such a way that the 
quadratic part of the transformed Hamiltonian function 

H (2 \z) = H (1 \M- 1 z) (143) 

assumes the form 

^ 2 W) = £f (Pl + ql). (144) 

k=l 
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For the first two steps in the sequence (|133p . we actually did not give explicit 
expressions for the generating functions W\ and Wi. For conceptual reasons (and to 
justify the notation) it is worth mentioning that such expression can be determined 
(see [Sj). The next steps in (|133p though rely on the explicit computation of the 
generating functions W n with n > 3. To this end it is convenient to introduce the 
adjoint operator associated with a phase space function A: 

ad A :B^ ad A B = {A, B} . (145) 

The transformation fj 134j) then leads to a transformation of the Hamilton function 
^•(n-i) £ Q jj(n) n > 3 which ~ m terms of the adjoint operator reads 

00 ^ 

H {n) = Y.^w n \ k . (146) 

k=0 

In terms of the Taylor expansion defined in Eqs. (|139p - (|14ip the transformation 
introduced by Eq. (|146p reads 

E y^w n ] k H^ k H 2) , (147) 

k=0 

where |_ - J gives the integer part of a number, i.e., the 'floor'-function. 

Using Eq. (|147p one finds that the transformation defined by (|146p satisfies the 
following important properties for n > 3. Firstly, at step n, n > 3, the terms of order 
less than n in the power series of the Hamiltonian are unchanged, i.e. 

W = #("-1) , for s < n , 
so that, in particular, H% = . Defining 

P = ad {2) ={Hf\-} (149) 

we get for the term of order n, 

tfW = - VW n . (150) 

This is the so-called homological equation which will determine the generating func- 

(n) ■ (n) 

tions W n for n > 3 from requiring T>Hn = 0, or equivalently Hn to be in the 
kernel of the restriction of T> to W n . In view of (|150p this condition yields 

R (n-1) _ e Ker p| w?l _ ( 151 ) 

Section 3.4.1 of Ref. [44J provides the explicit procedure of finding the solution of 
Eq. (|15ip . In the generic situation where the linear frequencies u)\ , . . . , ojf in (|144|) 
are rationally independent, i.e. m\UJi + . . . + rnjujf = implies mi = . . . = m/ = 

(n) 

for all integers mi, . . . , mf, it follows that for odd n, Hn = 0, and for even n, 

H n n ^ e span [ I^I^I^ 3 . . . f/ : \a\ = n/2} , (152) 



(148) 
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where I k = (q\+p\)/2, with k = 1,... ,/. 

Applying the transformation ()146p , with the generating function defined by (j!50|) , 
for n = 3, ... ,no, and truncating the resulting power series at order no one arrives 
at the Hamiltonian H^ 1 ^ corresponding to the n^ 1 order normal form (NF) of the 
Hamiltonian H: 

n 

H^\z)=E + ^2h^(z). (153) 

s=2 

The normalized Hamiltonian H^t is an Uq 1 order approximation of the original 
Hamiltonian H obtained from expressing H in terms of the normal form coordinates 
znf which in turn are obtained from the symplectic transformation of the original 
coordinates z = (qi,Pi) 

z NF = <p(z) = (<f>w„ 0W„ o -i • • • 4>W 2 <?Vi)( z ) • (1 54 ) 

This is why one can use instead of H to analyze the dynamics in the neigh- 

borhood of an equilibrium. 

The procedure to compute JSjjjp and the corresponding coordinate transforma- 
tion is algebraic in nature, and can be implemented on a computer. A computer 
program is freely available from [46J. 
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